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Abstract 

We present a model of optical communication system for high-bit- 
rate data transmission in the nonreturn-to-zero (NRZ) format over 
transoceanic distance. The system operates in a small group velocity 
dispersion regime, and the model equation is given by the Whitham 
equations describing the slow modulation of multi-phase wavetrains 
of the (defocusing) nonlinear Schrodinger (NLS) equation. The model 
equation is of hyperbolic type, and certain initial NRZ pulse with 
phase modulation develops a shock. We then show how one can ob- 
tain a global solution by choosing an appropriate Riemann surface on 
which the Whitham equation is defined. The present analysis may be 
interpreted as an alternative to the method of inverse scattering trans- 
formation for the NLS solitons. We also discuss wavelength-division- 
multiplexing (WDM) in the NRZ format by using the Whitham equa- 
tion for a coupled NLS equation, and show that there exists a hydro- 
dynamic-type instability between channels. 
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1 



1 Introduction 



Recently there has been a great deal of research works on designing a long- 
distance and high-speed optical communication system for the next genera- 
tion of lOOGbit/s over 9000km (transpacific distance). In the fall of 1996, 
the undersea optical cable called TPC-5 between Japan and US has been in- 
stalled, and the system operates at lOGbit/s in a non-return-to-zero (NRZ) 
format. Also in 1998, a total of lOOGbit/s (e.g. 5Gbit/s x 20Channels) 
wavelength-division-multiplexed (WDM)-NRZ system will be set up as TPC- 
6 optical cable. In such a system, noise accumulation and the combined ac- 
tion of fiber dispersion in the group velocity of pulse and fiber nonlinearity 
called the Kerr effect are the main system limiting factors. The analysis 
of NRZ signal propagation was hampered so far by necessity of lengthy nu- 
merical simulation and expensive experiments. Recently, we proposed in 
a hydrodynamic model to describe nonlinear effects in the NRZ pulse 
propagation, which was derived as a weak dispersion limit of the nonlinear 
Schrodinger (NLS) equation. Based on this model, we also discussed meth- 
ods to control the pulse broadening by means of initial phase modulation 
(pre-chirping) |jl5[ and nonlinear gain . 



In this paper, we give a mathematical formulation of the model in the 
framework of the Whitham averaging method which describes the slow 
modulation of the amplitude and phase of quasi-periodic solution of the NLS 
equation. We start in Section 2 to provide a necessary information of the 
problem, and derive a shallow water wave equation as a simplest approxima- 
tion of the model. In Section 3, we give a mathemtical background of the 
Whitham equation as describing dynamics of the genus g Riemann surface 
determined by the spectral curve of the Lax operator for the NLS equation. 
Here the genus represents the number of periods in the quasi-periodic solu- 
tion, and the Whitham equations are given by first order quasi-linear system 
of partial differential equations. In particular, we identify the shallow water 
wave equation as the genus zero Whitham equation. We also note that the 
spectral curves corresponding to the Whitham equation for the NLS equation 
(NLS- Whitham equation) are hyper-elliptic. In Section 4, we first show the 
hyperbolicity of the NLS- Whitham equations which enables us to classify the 
initial data for the global existence of the solution. We then decribe the effect 
of pre-chirping for reducing the broadening of NRZ pulse, and as a result we 
find that the Riemann surface becomes singular (formation of a shock) , and 
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its genus changes from zero to either one or two depending on the strength 
of the chirp. 

Since the NRZ system operates at a small group velocity (i.e. second 
order) dispersion regime, the third order dispersion becomes important, es- 
pecially for pulses with shorter width (in a higher bit rate system). We 
discuss this effect in Section 5, and show that the third order dispersion 
leads to asymmetric distortion of the NRZ pulse and results a formation of 
shock-like structure. In Section 6, we extend the NLS-Whitham model to de- 
scribe pulse propagation in a WDM-NRZ system. The extension is obtained 
by taking a weak dispersion limit of a coupled NLS equation, in which the 
interaction between channels due to the Kerr nonlinearity is expressed by 
the cross phase modulation. We then find a hydrodynamic-type instability 
between channels, which corresponds to loosing the hyperbolicity of the hy- 
drodynamic model obtained in the dispersionless limit of this coupled NLS 
equations. The WDM model is further studied by using an integrable model, 
the Manakov equation or the vector NLS equation. The Whitham equa- 
tion of genus zero in this model is the well-known Benney equation which 
describes shallow water waves in a stratified fluids [@, and the number 
of channels in WDM corresponds to that of the layers. We then present how 
the nonzero genus Whitham equations appear in the Manakov equation. The 
main feature of this model is that the spectral curve is non-hyperelliptic, and 
the corresponding generic Riemann surface has a topology of a sphere with 
3h handles with some positive integer h. 



2 Preliminary 

For a long-distance and high-bit rate optical communication system, both 
dispersion and nonlinearity of the optical fiber are important and the NLS 
equation has been used as a model equation expressing such effects [^j. 

Here q is the complex envelope function of the electric field in the fiber, Z 
the propagation distance, and T the retarded time. The coefficients P2 and 
z/(> 0) represent the group velocity dispersion (GVD) and the Kerr (nonlin- 
ear) effect of the fiber. Here we further normalize q so that we can set z/ = 1. 
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(a) Bright Soliton 

(b) Dark Soliton 
(c) NRZ 

' — \ I \ I — \ I \ 

0010110010111100 

Figure 1: Various types of a 16bits coded signal (0010110010111100). 

It is well known that if /32 > (focusing or modulationally unstable case), the 
NLS equation has so-called bright solitons with sech-shape, and if /?2 < (de- 
forcusing or modulationally stable case), so-called dark ones with tanh-shape. 
These are the stationary solutions for Eq. (|2.1|) , and Fig.[l| shows examples of 
pulse sequences in |gp as a 16bits-coded signal (0010110010111100) using 
those pulses. In Fig.|l|, we also show the coded signal using the NRZ pulses. 
Because of the modulational stability of constant amplitude, for example, in 
the (011110) coding part, one has to operate the NRZ system at the deforcus- 
ing regime, (32 < 0. The NRZ pulse is however non-stationary, and in order to 
reduce the pulse distortion the system is designed in a small dispersion and 
low power regime. In an ideal case, the best performance can be obtained 
in the zero dispersion and linear limits. However, the real system always 
includes a noise, and the signal pulse power should be kept more than the 
noise level. Because of this limit in signal-noise ratio (SNR), the nonlinearity 
becomes an essential effect for a long distance transmission problem, and the 
zero dispersion causes a resonant interaction process between noise and sig- 
nal, which is called the four wave mixing (FWM), and leads to a distortion 
of the pulse. The main purpose of this paper is to develop a mahematical 
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theory to describe the pulse distortion in a weak dispersion hmit of the NLS 
equation. Here we do not discuss, however, the distortion due to the FWM 
with noise. 

Let us first recall that the NLS equation (|2.1| ) has an infinite number of 
conservation laws, for example, 

^ V|gp4Arg(g)V (2.2) 




dT \ dT 



^1 df^^^(^n = df[-Tdf^^^ 



dq 



dT 



+ \\^\"\- (2-3) 



Here we have set [32 = — The quantities |gp and Arg(g):= (i/2) ln(g*/g) 
represent the local intensity and phase of the electric field. We then introduce 
p := \q\ and a := /3Arg(g), and rewrite q as 



g(T, Z) = ^p(T, Z) exp U-a(T, Z) 1 . (2.4) 



Substituting Eq.( p.4|) to Eqs.( p.2|) and (|2.3| ), and rescaling the variables Z — *• 
Z/P, we have 

i^P = 4F (P«) > (2-5) 



dZ"^ dT 



where u represents the chirp defined hj u := da/dT. In a small dispersion 
limit (3 0, if we assume that both p and m are smooth, then Eqs.( p^.5[ ) and 
( |2.6| ) for p 7^ can be approximated by a 2 x 2 quasi-linear system, 

d r p\ _ f u p\ d f p-^ (^2 7) 



dZ \uj VI uj dT \u 

A detail on the weak dispersion limit of the NLS equation has been studied 
in [0, where the existence theorem of the limit is given based on the Lax- 
Levermore theory. When p > in ( |2.7|) , the eigenvalues of the coefficient 
matrix are real, u + ^p and u — ^/p, and the system is totally hyperbolic. 
The system ( p.7| ) is known as the shallow water wave equation, and has been 
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intensively discussed (for example, see [^). It is then interesting to note 
that the distortion of the NRZ pulse may be understood as a deformation of 
the water surface. 

For a demonstration of the NRZ pulse propagation, we consider the initial 
value problem of Eg. ( p.7|) . As a simple example of the initial NRZ pulse, we 
take a square pulse having constant phase (zero-chirp). 



P(T, 0) 



/ Po, 


for 


|r| 




I 0' 


for 


T 





u{T, 0) = 0, for VT. (2.8) 



This is called "Dam-breaking problem", since p and —u represent the depth 
and the velocity of water which rests on the spatial region \T\ < Tq at the 
time Z = 0. We then expect to see a leakage of the water from the edges, 
and in fact we obtain: 

Proposition 2.1 The system (^. 7\ ) with the initial data has the fol- 

lowing global solution up to the distance Z = To/^/po-' ^or T > 0, 



piT, Z) 
piT, Z) 
<T,Z) 



u{T, Z) = 
1 



for T > To + 2 Vp^Z, 



mzn 



Po, 
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-max 



0, 



2v^- 
Po + 



T-T 



Z 

T-To 



for <T <To + 2^Z, 
for 0<T <To + 2y%Z, 



For T < 0, we have p(T, Z) = p{—T, Z) and u{T, Z) = —u{—T, Z). 

Proof. As in the standard analysis of the method of characteristics for 
the quasi- linear hyperbolic system ( |2.7| ), we first rewrite the system in the 
Riemann invariant form p5| 



dZ VA2 



1 /3Ai + A2 



4 V 







d /Ai 







Ai + 3A2/ dT \\ 



where the Riemann invariants Ai and A2 are given by 

Xi = u- 2^, A2 = M + 2y/p. 



(2.9) 



(2.10) 
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Then the initial data for the Riemann invariants are given by 



Ai(T,0) 



A2(T,0) 



for T < -To 
-27p^, for T>-ro ' 

(2.11) 

for T < To 
-2^/p^, for T > To ■ 



Here we extend the values of Ai and A2 for the region |T| > Tq. One should 
note from Eq.( p.6|) that u is defined only the region where p 7^ 0. Then 
Eq. (|2.9|) becomes a single equation for A2 near T = Tq where Ai = 

i(-2VA; + 3A2)^, (2.12) 



from which we have the solution expressing a rarefaction (centered simple) 
wave starting from T = Tq; for T > 0, 



2Vp^, for < T < To - 

A2(T,Z) = <| ^^_1^^_Z£, for To-^Z<T<To + 2^Z, 

-2^, for T > To + 2^Z. 

(2.13) 

This solution is valid up to the distance Z = Tq/^J])^ where the centered 
simple waves from T = To and — Tq collide at T = 0, and the solution should 
be modified. To find Ai(T, Z), note that the point symmetry in the graphs 
(Afc,T), i.e. A2(T) = — Ai(— T). Then from Ai and A2 we have the solutions 

P=^(Al-A2)^ « = l(Ai + A2). (2.14) 

Noting the symmetry p{—T,Z) = p(T,Z) and u{—T,Z) = —u(T,Z), we 
complete the proof. □ 

The evolution of the Riemann invariants Ai and A2 is shown in Fig.^. In 
Fig.|^, we plot the pulse shape p and the chirp u obtained in Proposition 



The thin dotted curve in Fig.y shows the numerical result of the NLS equation 
with po = l,To = 10 and (32 = -p^ = -0.1 at Z = 10/3 ^ 3.16. Notice the 
good agreement with the analytical solution except some small oscillations 
on the top of p which disappear in the limit (32 —>■ 0. As we have predicted 
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Figure 2: Evolution of the Riemann invariants. Dashed curves are the initial 
data: Solid curves are the solution Afe(T, Z) for Z > 0. 




Figure 3: Deformation of NRZ pulse (Dam breaking problem). The thin 
dotted curve on p shows the numerical solution of the NLS equation. 
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Figure 4: Optical shock due to initial chirp. 

from the hydrodynamic analogy, the deformation of the pulse induces the 
generation of chirp (or water velocity) at the edges, and the global nature 
of the solution is understood as an expansion of the water (the rarefaction 
wave). The expansion is characterized by the characteristic velocity s at the 
edge, and for example at T = Tq it is given by s = — i/po for T < Tq and 
s = I-s/pq for T > Tq. In order to reduce the expansion, it is natural to put 
initial chirp opposed to the chirp appearing in the edges. This is equivalent 
to give an initial velocity with a piston and to compress the water. Because 
of the quasi- linearity of the system (|2.71 ), we then expect a shock formation 
in the solution, and therefore Eq. ( |2.71) is no longer valid as an approximate 
model of the NLS equation. In fact, if we start with the initial data with the 
same p in ( |2.8| ) but nonzero chirp u, for example, with uq > 

the numerical solution of the NLS equation results high oscillations starting 
from the discontinuous point of chirp, T = 0, as seen in Fig.§. Here we 
have set po = 1,mo = 3/3 ~ 0.95 with P2 = —0.1, and the solution p{T,Z) 
is at Z = 5/3 ~ 1.58. This phenomena has been also noted in |^ as an 
optical shock appearing in somewhat similar situation. The main frequency 
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of the oscillations is of order 1/(3, and usual optical filter can remove those 
oscillations as in the sense of averaging. So what we want to describe here is 
the average behavior of the solution. For this purpose, in the following section 
we extend the model (p.7|) and show that the new model admits a global 
solution describing the average motion of the NLS equation (|2.1| ) for several 
step initial data having different values of chirp. This is the main purpose 
of the present paper, and it may be referred as a dispersive regularization 
problem of the shock singularity of hyperbolic system in a weak dispersion 
limit of the NLS equation. 

3 The NLS-Whitham equations 

We here develop a regularization method based on the integrability of the 
NLS equation. The key idea is to identify the shock singularity as a gener- 
ation of a quasi-periodic solution of the NLS equation which appears as the 
high oscillations as seen in Fig. 4. The quasi-periodic solution can be char- 
acterized by the Riemann surface defined by the spectral curve of the Lax 
operator in the inverse scattering transform (1ST). Then the genus of the 
Riemann surface changes through the singularity. Here the number of the 
genus corresponds to that of independent phases in the quasi-periodic solu- 
tion 1^, Then the main purpose here is to derive an equation to describe 
the dynamics of the Riemann surface corresponding to the modulation of the 
quasi-periodic solution of the NLS equation. The resulting equation is called 
the NLS-Whitham equation. 

Let us first recall the 1ST scheme of the NLS equation and reformulate the 
quasi-linear system (|2.7|) as the simplest approximation of the NLS equation 
in a weak dispersion limit, that is, the genus zero NLS-Whitham equation. 
Then we discuss the general case with arbitrary genus. Here we somewhat 
use a different formulation of the 1ST method, which is related to Sato's 
formulation of the KP hierarchy [^. This is quite useful not only for the 
present purpose but also for an integrable extension of coupled NLS equation 
used as a model of WDM system discussed in Section 6. The Lax operator 
L of the NLS equation in this formulation is given by a pseudo-differential 
operator with D := —ifid/dT, 

L = D + ^p{D-u)-^^p 
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D + pD-^ +[pu- ]-Dp ) + . . . . (3.1) 



2 

This can be derived from the standard Lax operator in a 2 x 2 matrix form 
[ p^ , and the operation with is given by a generahzed Leibnitz rule with 
DD-^ = D-^D = 1, 

^'/■ = f:(f)pV)^'-'- , for anyfeZ. (3.2) 

Then the 1ST scheme (the Lax pair) for the NLS equation is given by the 
eigenvalue problem Lcj) = Xcf) and the evolution equation, 

-^11 = 520, B, = \{L%_,:=\d' + p. (3.3) 

Here the symbol (i^^)>o denotes the differential (i.e. polynomial in D) part 
of the pseudo-differential operator L^. 

As a genelarization of this scheme, the integrable hierarchy of the NLS 
equation can be formulated as 

-z^ = 5„0, 5„ = -(L")>o, (3.4) 

where is the parameter describing the evolution generated by i?„ for 
n = 1,2, In particular, Zi represents the translation in T, that is, 
0(r, Zi, Z2, ■ ■ ■) = 4>(T + pZi, Z2, ■ ■ ■) so that we identify pZi as T, and Z3 
represents the parameter for the complex modified KdV equation which will 
play an important rule in Section 5 where we discuss a higher order correc- 
tion to the NLS equation. The compatibility conditions for the 0-equations, 
L(f) = \(f) and (|3.4|) , then give the Lax equations 

dL 

- I—- = [5„, L] := BnL - LBn, (3.5) 

OZn 



from which we recover Eqs.( |2.5| ) and ( |2.6| ) for n = 2. It should be noted here 
that all the flows in (|3.4| ) are compatible, and this is a direct consequence 
of the definition of i?„ [^. Namely we have: 

Lemma 3.1 The differential operators i?„ in (\3.4\ ) satisfy the following re- 
lation called the Zakharov-Shabat equation, 

-^|^+^|^ + [S«^,5n] = 0, (3.6) 

which gives the compatibility conditions among the flows in ( ^.d{ ). 
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Proof. Consider the equation, 

Then sphtting the operator L™ into the parts having nonnegative and neg- 
ative powers of D, i.e. = + with B^ := (L'^)<o, the right hand 
side can be written in the form, 

Projecting the equation on to the nonnegative (differential) part, we obtain 
the result. □ 

Now we take the limit /? — which corresponds to a quasi-classical limit 
of the NLS equation, and assume to be of the WKB form, 

= exp {^-^S^ (3-7) 

where S is called the action. Then, noting D'^^/^ {dS/dTy as /3 — for 
any integer z/, the eigenvalue problem L0 = A0 in the limit gives an algebraic 
equation for the spectral A, 

A = P + -^, (3.8) 
F — u 

where P is the momentum defined by 

P := ^ = lim f -D(j)] . (3.9) 

The time evolution ( |3.3| ) with the rescaling Z Z/(3 gives 

dP dQ2 



dZ dT 



(3.10) 



where Qn is defined by lim^_^o(-Bn0/0), which is just the polynomial part of 
A"/n in P. We denote this as Qn = [A"]>o/n, and for example Qi = [\]>o = 
P^Q^ = [A2]>o/2 = + P- Equation (|3.10|) then provides the system 
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( |2.7| ). In particular, we note that the momentum P gives the (polynomial) 
conserved densities P„ = Pn{p,u) of the system (|2.7| ), i.e 



oo 



P^.-Y.^^^-i-^-: (3.n, 



This is also true for the original NLS equation in the Lax form ( p.5|) . Namely, 
first express D in the Laulent series of L from ( p.lD , 



D = L-J2 DnL-"" = L- pL-' -(pu-^Dp)L~^~---, (3.12) 



n=l 



where D„ are the differential polynomials of p and u, which we denote Dn = 
Dn[p, u : P]. Also from Eq.( p.3| ) we have 

which shows that the D„ are the conserved densities of the NLS equation. 
The dispersionless limit of Dn then gives P„, 

P^ip,u) = \imD^[p,u:/3]. (3.14) 

Note here that all the derivatives of p and u in Dn vanish in this limit. This is 
valid if the functions p and u are both smooth in T. However, as seen in Fig.^ 
these functions develop high oscillations for certain initial data. Then one 
should extend the system ( |2.7| ) to the case including those oscillations. This 
is indeed the main purpose of this paper, and the extended system is called 
the Whitham equation describing the averaged behaviour of the solution. 

We also note from Eq.( p.4|) that the hierarchy can be reduced to the 
following form with Z„ —>■ Zn/ (3 and Zi = T, 

dQn dQni ^2 ^^-^ 

dZni dZn 

which give the compatibility conditions of the dispersionless limit of 

Qm. (3.16) 



dZ„i 
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Note that Eg. ( p. 161) with P = dS/dT gives the Hamilton- Jacobi equation 
with the hamihonian Qm = Qm{T,dS/dT, Zn)- This quasi-classical limit 
has been discussed in the framework of the dispersionless KP hierarchy (see 
for example |l|, 0). 



From this formulation, we now give a geometric interpretation of the 
system ( |2.7|) . First note that the algebraic equation ( |3.8|) determines a two- 
sheeted Riemann surface of A with branch points Xi = u — 2y/p and A2 = 
u + 2y/]). Then after compactifying each sheet and gluing them together 
along the branch cut, we have a Riemman surface of genus g = which is 
defined by the curve y'^ = Ro{\) = (A — Ai)(A — A2). Namely Eq.( |3.8| ) can 
be written in the form [P — (A + m)/2]^ = and the roots are given by 

1 



P=-[^\ + u±^R,{X)j. (3.17) 

The branch points Ai and A2 are the Riemann invariants of the hyperbolic 
system and the values Pi = P{\i) and P2 = -P(A2) give the characteristic 
velocities which are also given by the roots of the equation dX/dP = 0. 
The system ( p.7|) then describes a (slow) modulation of the spectral which is 
invariant for the NLS equation. In Fig.^, we plot the curve ( |3.8| ) and those 
points. Note that the spectral of L is given by Spec(L) = (—00, Ai] U [A2, 00). 
The g = Riemann surface has a topology of a sphere as shown in Fig.^ 
Thus the hyperbolic system (|2.7|) can be considered as an approximate system 
of the NLS equation on the g = Riemann surface. Let us first demonstrate 
that the system (|2.7|) or (|3.1CI| ) can be interpreted as a dynamics of the 
Riemann surface of genus 0. On the g = Riemann surface, we define the 
following meromorphic (Abelian) differentials ui and uj2 of the second kind, 

A2 




c^i = - ( 1 + ^^ ]dX^dX + 0[ — \ , (3.18) 



2 V v-Ro / 



AdA + 0(^), (3.19) 



where ai = Xi + A2(= 2m) and o"2 = A1A2. Note from Eq.( p.l7|) that these 



differentials are given by cui = dQi = dP and cu2 = dQ2 = PdP, that is, Qi 
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Figure 5: The g — algebraic curve (3.8). 




and Q2 are the Abelian integrals. In the case of (7 = 0, these differentials are 
trivial (image of the derivation), and it corresponds to that the fundamental 
group for = surface consists of only identity, i.e. for any cycle a on the 
surface, 

(fuj^ = 0. (3.20) 

Thus we have: 

Proposition 3.1 The quasi-linear hyperbolic system (^. !( ) can he written by 

duji duj2 



dZ dT 



(3.21) 



where uji and lo^ ore the Abelian differentials ^3.18\) and ( \3.19i ) defined on 
the g = Riemann surface. 

Equation ( p. 211) then defines the {g = 0) Whitham equation [|, The 
differentials Ui and uj2 in Eq.( p.21j ) are also given by the limits. 



uj-i = dQ-i = d lim ■ 



UJ2 = dQ2 = dihm --B2C, 

\ /3^0 (p 



(3.22) 
(3.23) 



Thus Eq.( ^.21D is just the compatibility condition of the T- and Z-flows, that 
is, d'^\n(t)/dTdZ = liKp /dZdT. Calculating the residue of (IQiP at each 
branch point A^, we obtain the Riemann invariant form of ( p.21|) . 



Sk{M, A2) 



A; = 1,2, 



dZ " 

where the characteristic velocity Sk{\i, X2) is given by 



Sk 



A 



(3.24) 



(3.25) 



A=Afc 



Equation (|3.24| ) is of course the same as ( p.9|) . Then we can rephrase the 
Proposition 2.1 as: 

Theorem 3.1 The g = Whitham equation l \3. 24\ ) with the initial data 
l\2.1[\) has a unique global solution (rarefaction wave). 
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Figure 7: The g = I algebraic curve = Rg{X) in (|3.26| ). 



We now present the NLS-Whitham equations for arbitrary genus. The 
two sheeted Riemann surface of genus g is defined by the algebraic (hyper- 
elliptic) curve y'^ = Rg{X) with 

2g+2 

R9W= Ei^-^k), (3.26) 

k=l 

where the branch points Xk are real and are assumed to satisfy Ai < A2 < 
. . . < X2g+2- The Riemann surface of genus g has a topology of a sphere with 
g handles. In Fig.|^, we plot the curve of g = 1. Comparing the case with 
g = 0, we see that the closed curve in the middle of the graph indicates an 
opening of the genus. 

The NLS-Whitham equation on the genus g Riemann surface is derived 
from the NLS equation as the slow modulation of quasi-periodic wavetrain 
where the spectral of the Lax operator L is given by, 

Spec(L) = R\{U(A2.+i, A2.+2)}. (3.27) 

Here we reformulate the result in as an extension of the g = NLS- 
Whitham equation ( p.21| ). In the case of arbitrary genus g, the equivalent of 
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Figure 8: The 2-sheeted Riemann surface with branch cuts and canonical 
cycles. 



the Abelian differentials ui and 002 defined before are given by 

1 / - + aiA^-i + ■ ■ ■ + . ,^ 
co-i = - IH 1= 1 d\ 



0J2 



dX + 



A2 



as A — s> 00, 



(3.28) 



1 
2 



A + 



A^^+2 - iaiA^+i + (i(T2 - l(rf)\3 + + ■ ■ ■ + Pg 



dX 



'dX 



XdX "'"^(^j' A— >oo, 



where ai, Oi (and (T3) are defined by 



(3.29) 



i<i i<j<k 



(3.30) 



The coefficients a„ and /?„ in ( |3.28|) and ( p.29| ) are determined by the nor- 
malizations on Ui and UJ2, 



U)i = 
ai J a. 



uj2 = 0, for i = 



(3.31) 



with the canonical aj-cycle over the region [A2j+i,A2i]. Fig.reff3:l shows the 
spectral of L and the canonical cycles. Namely, we have the following 2g 



18 



linear system of equations for at and (3i hi i = 1, ■ ■ ■ , g, 



If^' - -ajf + ajf-' + . . . + = 0, 

/I 1 \ (3-32) 
^iir' + (2^2 - -a!) If + f3jr + ■ ■ ■ + = 0, 



7-9+2 

2 



where the integrals /" are given by 

1 rMi+2 \^ 



d\. (3.33) 



Then the NLS-Whitham equation of genus g is defined by the same form 
as Eq.( p.21 ). As shown in ||^, 18, the Whitham equation describes the 



conservation laws averaged over the high oscillations appearing as a dispersive 
shock. Then one can introduce the fast and slow variables to describe the 
oscillation and the averaged behavior on the solution of the NLS equation 
in the weak dispersion limit ^. The average is then taken over the fast 
variable so that the averaged equation describes the solution behaviour in 
the slow scale. Thus the average is given by 

if) = {f){Tsiow) ■= hm— /" f{Tsiow,Tfast)dTfasu (3.34) 

where Tgiow and Tfast are the slow and fast variables, and we use T as the 
slow variable in the NLS-Whitham equation. Note that the rescaling Z — >■ 
Z//3 used to derive the Whitham equation (|3.10|) defines the slow variable 



Z = Zsiowi i-e. Zsiow = f^Zfast- Then the Abelian differentials ooi and uj2 in 
the Whitham equation (|3.21|) are defined by the averages 0], 

^1 = < - In 0) ~ (^1 + E ^ j dX, (3.35) 

u^2 = d{- A In 0) ~ (^A + £ ^) dX. (3.36) 

The Whitham equation leads to the conservation laws, 

dFn dGn 



dZ dT 



(3.37) 
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in which the cases for n = 1 and n = 2 correspond to the averages of ( |2.2| ) 
and ( p.3|) . Here F„ and Gn are the averaged conserved densities and fluxies, 
and we have for example 

^1 = (P) = ~\^2 + ^ctI + ^ai, (3.38) 

I 1 1.1 

F2 = Gi = (pu) = -(73 - gf^lf^2 + —(yf + -Pi 

II 1 / N 

= gC^s - g^Tias + ^(Ti + -aiai + -as, (3.39) 

^ / 2 1 2\ 1 1 12 5 2 

Gs = + -p ; = --a, + -a,a, + -a^ - -a,a, 

+ + l^iPi + \P2. (3.40) 

One should note that the quantities F„ are given by the averages of Dn in 
( |3.12|) , i.e. Fn = (Dn). In the particular case of 5^ = 0, we have ai = Pi = 0, 
and Fn become P„ in ( |3.11| ), which means there is no oscillations in these 
quantities. 



4 Control of the NRZ pulses 

We now consider the initial value problem of the NLS-Whitham equation 
with the initial pulse having a nonzero chirp in a purpose of reducing the 
pulse broadening due to nonlinearity and the dispersion. The main objective 
is to determine the average behavior of the solution of the NLS equation ( |2.1] ) 
in the weak dispersion limit. In a practical situation, this describes the NRZ 
pulse behaviour through an optical filter which removes the high oscillations 
(shock) due to the initial chirp. 



4.1 The hyperbolic structure of the NLS-Whitham equa- 
tion 

Let us first discuss the hyperbolic structure of the NLS-Whitham equation, 
which enables us to analyze in detail the initial value problem. As in the case 
of (? = 0, we have the Riemann invariant form by evaluating the residue at 
A = Afc in the NLS-Whitham equation (|3.21| ) with ui and 002 given by ( p.2^ ) 
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and ( p:29D , 



Sfc(Ai,---,A2g+2)^::T^, tor k 



2g + 2. 



(4.1) 



where the characteristic velocities Sk is defined by the ratio of the polynomials 
in cui and UJ2, i-e. 



Af+i 



|(TiA9 + aiA9-i + ■ ■ ■ + ttg 



(4.2) 



A— At 



Then following the ref.|T^, we obtain: 
Lemma 4.1 The characteristic velocities satisfy 

dsk 



ii) 



> 0, 



for Wk, 
if Xj > Afc, 



(4.3) 
(4.4) 



dXk 

Sj > Sky 

Proof. Let us call the polynomials in Eq.( |4.2| ) P(A) for the denominator 
and Q{X) for the numerator. Then from the normalizations ( |3.31| ), we have 

• -P(A) has exactly one zero in the gap [A2j+i, A2i+2] for z = 0, 1, 

• Q{X) has at least one zero in each gap. 

Q{X) 



.9, 



s{X) 



P(A) 



X + 



as A 



00. 



From these proparties, the following function can be expressed as the curve 
in Fig.|; 



y,(A) := six) - s, - 
It is then easy to see the relations 



^ (g(A) - S.PW) . 

and (|0). □ 



(4.5) 



From this Lemma, we have the following useful Corollary which provides 
a way to regularize the problem for the global existence of the solution: 

Corollary 4.1 // all the initial Xk are monotonically decreasing and sepa- 
rated in the sense that for VT G R and VA; G {1, 2, ■ ■ ■ ,2g + 2} 



mm[Afc+i(T,0)] > A,.(T,0) > max[Afc-_i(T, 0)], 
then the NLS-Whitham equation has a global solution. 



(4.6) 
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V2k.l(?l) 




4.2 Regularization of shock singularity 

We now study the initial value problem of the NLS-Whitham equation with 
the initial data shown in Fig.p!0|, i.e. 

p(T 0) = I !^! < 

' ' \ 0, for |T| > To ' 

(4.7) 

r -Mo, for -To < T < 
m(T, 0) = < uo, for < T < To . 
[ 0, otherwise 

Here we vary the value of initial chirp uo, and study its effect for a purpose 
of reducing the NRZ pulse broadening. As an example of real situation, we 



show in Fig.|TT] a numerical result showing the effect of initial chirp for a 
16bit coded data (0010110010111100). Here the initial phase modulation is 
set to be periodic with the period given by the width of one-bit pulse (bit- 
synchronous modulation), and the output is smoothed out by using a filter 
(averaged over high oscillations). As seen in the Figure, the overlapping 
of pulses in a) with no initial phase is suppressed by the periodic initial 
phase modulation. The period is seen for example in (00111100) part as the 
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-To 



To 



Uo 



-To 



T 







To 



-Uo 



Figure 10: The initial phase modulation for one-bit NRZ pulse. 



elevation and depression of the pulse level. It looks that NRZ pulse with 
an initial phase modulation tends to deform into an RZ pulse which may 
have a better property in signal processing in a network system of optical 
communication. 

Before we start to analyze this solution behaviour in detail, let us recall 
that the value of the u corresponds to the initial velocity of the water, and the 
positive (negative) value of u gives an compression (expansion). So we see 
a shock formation for the positive u and a rarefaction wave for the negative 
u, which explains the pulse deformation shown in Fig.|ll]. We also note from 
Eq.( |3.1| ) that the characteristic speed for zero chirp is given by 2y/p, so that 
we expect to see a different behavior of the solution for each uq value, {i) 
uq > (a) < Uo < 2y/po, {iii) > mq > — 2y/po, and {iv) uq < —l^fp^. 

This situation is quite similar to the case of Toda lattice equation discussed 
in 1^. We then show that by choosing an appropriate number of genus for a 
given initial data and by solving the corresponding NLS-Whitham equation, 
one can obtain a global solution. In the form (|4.7| ), we here consider the case 
with To = oo, since the edges of the pulse with constant chirp do not give 
any essential change from the case with zero chirp except the velocity of the 
expansion. 

We now consider each of the four cases separately: 
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Figure 11: Effect of initial phase modulation for NRZ pulse propagation, a) 
Without, and b) with the initial modulation. 



(i) The case where Uq > 



In this case, the g = Riemann invariants Ai and A2 at Z = are given 
by increasing functions, i.e. 



Ai(T,0) 



A2(T,0) 



-Mo - for 
Mo - 2y^, for 



-no + for 
Uq + 2^, for 



T < 

T > ' 

T < 

T > ■ 



(4. 



Then we see from Corollary [4.1| that the g = NLS-Whitham equation de- 
velops a shock. So we have to use the nonzero genus NLS-Whitham equation 
for the regular solution. First we note: 

Lemma 4.2 The initial data can be identified as the following g = 1 data, 
Ai = -Mo - 2v^, for VT, 




-uo + 2^, for T < 0, 
-Mo - 2^, for T > 0. 
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(4.9) 




Mo + 2^, for T < 0, 
for T>0, 



A4 = uo + 2^, for \fT, 



Proof. To prove the Lemma, we need to show the values (p) and (pu) 
in (|3.38|) and (|3.39| ) for (7 = 1 become (p) = p and (pw) = pw at the point 
Z = 0, that is, the initial data is a degenerate g = 1 data. The quantities ai 
and /5i in ( |3.38| ) and ( |3.39D are determined by the normalization (|3.31|) , that 
is, from ( p.32| ) 



with 



- -aJl + aJl = 0, 



- ) il + A/? = 0, 



(4.10) 



Tk 1 

A = o 



A4 



2 7a3 /_(A-Ai)(A-A2)(A-A3)(A-A4) 



--d\. 



Substituting A = A3 + (A4 — A3) sin^ 6' as in [Q, we have a convenient form of 

Tk 

-'1 5 



If 



7r/2 



[A3 + (A4-A3)sin2 0]'= 



[A3 - Ai + (A4 - A3) sin^ 6] [As - A2 + (A4 - A3) sin^ 6] 



--de. 



Now let us compute (p) and {pu). We first note that the Riemann invari- 
ants has the following symmetry in the T-coordinate: 



Afe(-r) = -A5-;l-(T), for fc = l,---,4. 



(4.11) 



so that we have 



(Tk{-T) 

m-T) 



-l)Vfc(T), for k = l,-- 
-iriUT), for n = l,2. 



(4.12) 



25 



Thus we consider only the case with T < 0. We then have = X'^Ii and 
ai = — mo(mo + '^y/Po),Pi = 2po(Mo + S^ypo), indicating that we do have 
g = 1 initial data. Using these values, it is straightforward to show the result 
(p) = po and {pu) = — poMo for T < 0. □ 

The point here is to recognize that the initial data is degenerate and can also 
be identified as a non-zero genus data in which the Riemann invariants are 
decreasing functions of T. Then we have: 

Theorem 4.1 For uq > 2ypo, the g = I NLS-Whitham equation with the 



initial data ( U-^ ) ho,s a global solution. 



Proof. First note that =constant remains for \/Z. So we only need to 
determine the evolutions for A2 and A3. In order to compute these values, we 
first regularize them at T = 0: From the relation (|4.4| ) in Corollary [4.1| , we 
have S2(T) < s^iT) for VT, in particular, 52(0) < S3 (0) and s^(0) < S3'(0), 
where s^(0) := lim^-^Oit Sk{e). So if we impose 

A2(T) = lim Ai(T-e), and A3(T) = lim A4(T - e), (4.13) 

then we have sjj" > S2 , and therefore from ( [4.3| ) in Corollary ^?1| we obtain 
the global solutions for A2 and A3. We confirm this by caluculating explicitely 
the evolution of the Riemann invariants. The formula of Sk '.= s(A)|A=Afe for 
the g = I case is given by 

A"' — 2^iX + ai 



Using ( [4.10| ), this can be written by 



1 A'- J + (l^j-|^J)(A-it) 

^ A^--^-20-i(A 



^1 



We now compute S3 as well as S3", S2 , and S2". However, from the sym- 
metry discussed in the proof of Lemma we have S3 = — s^, S3 = — S2 , 
and also S3 > S3 from (|4.3| ) in Corollary So we need to compute only 
and s^, which can be evaluated for the region of T < 0. In the computation 
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of S3 , we have Ai = —A3 = — A4 = —Uq — 2y/po and A2 = —uq + l^fpo- Then 
we find l\ = Ag/^, and obtain 



In the case of S3 , we have Ai = — A4 = —Uq — 2y/po and A2 = —A3 = 
-Uq + 2y/pQ. Then we obtain 



^ _ -»g - 6uo^ + 4po + 6y/p^{uo - 2^)Ji + 8^J2 
~ Uq- 2^q{1 - Ji) ^ ■ ' 

where Jn,n = 1,2, are given by 
J, 



\l + a^sm'e){l + a^cos^e)' ^0 ^(1 + sin^ ^)(1 + cos^ ^) 

where = 2y/po/uo < 1. One can check the relation, 

S3 > s+ > 0. (4.18) 

This proves that the g = 1 NLS-Whitham equation with the initial data 
admits a global solution. □ 



In Fig.|T2|, we plot the evolution of A2 and A3 together with the constants Ai 
and A4. We also plot in Fig.|l3| the shape of the pulse p(T, Z) at Z = 3/5 ~ 
0.95 obtained by the numerical simulation of the NLS equation ( p.l|) with 
Po = 1, uo = 9/3 ~ 2.85 and P2 = —0.1. We see here nearly steady oscillation 
in the center region |T| < s^Z, which indicates the g = 1 periodic solution 
of the NLS equation as predicted in Fig.^. 

(ii) The case where 2^fp^ > mq > 0: 

This case also gives a compression of the pulse, and we expect a shock 
formation in the system Eq.(2/?). So we again need to use the nonzero genus 
NLS-Whitham equation for the regularization. As in the previous case, we 
first identify the initial data ( [4.7|) as a non-zero genus data expressed as 
decreasing functions: 
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Figure 13: Deformation of the pulse p{T, Z) for uq > 
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Lemma 4.3 The initial data Ij ) with 2y/po > Uq > can be considered as 
a g = 2 initial data, 



Uq - 2^, for T < 0, 
-Mo - 'iJJ^, for T > 0. 



Ai = -Mo - 2^/p;^, for VT 

A2 = 

As = Mo - 2v^, for VT, 
A4 = -Mo + 2v/a^, for VT, 

A5 = 



(4.19) 



mq + 2^, for T < 0, 
-Mo + 2yp^, for T > 0, 



Ag = Mo + 2yp^, for VT, 

Proof. The reflective symmetry in the data mentioned in the previous 
case holds here as well, and we have 

Afc(-T) = -A7-fe(T), for A;=l,---,6, (4.20) 

which lead to 

a,(-T) = (-l)Vfc(T), for A; = l,---,6. (4.21) 

The integrals /" of (|3.33|) with g = 2 also satisfy 

J|^(-T) = (-l)"Jf (T), for 2 = 1,2, and n = l,2,---. (4.22) 

As in the previous case, we need to compute (p) and {pu) for only T < 0. 
From ( |3.32| ), we see for i = 0,1 

= ([3,^]--a^[2,^]), (4.23) 



[1,0] V' ' 2 

-1)^+1 r., , 1 _ /I 1 
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where 

• Tn jm 
-'1 -'1 



det 



jn jni 
-'2 -'2 



For T < we note that Ai = — A5 = —Xq and A2 = A3 = — A4. Then 
caluculating [n,m] leads to ai = — (Af + A2)/2,/3i = (A2 — Ai)^/8, and using 
0-1 = (A2 - Ai), (72 = -{Xl + Xl + A1A2), erg = (Ai - X2){Xl + Xl), we verify 

(p) = ^(Ai + A2)' = Po, (4.25) 
{pu) = ^(Ai + A2)'(Ai - A2) = -pouo. (4.26) 

This completes the proof. □ 

From this lemma, we obtain: 

Theorem 4.2 For 2y/po > ^^o > 0, the g = 2 NLS-Whitham equation with 
the initial data ( U-i^ ) admits the global solution. 

Proof. We need to compute the velocities (= — s^) and sf(= —32)- 
The regularization of the initial data is the same as the previous case, and 
we here impose 

A2(T) = lim A2(T - e), X^iT) = lim X,{T - e). (4.27) 



Then calculating the integrals given by the normalizations ( |3.31|) , we obtain 
the foUwing results. 



S5 = -^2 = [uo + VPo) 1 + 



(% + v^)V' (4.28) 



4 = -■^2 = -^uo + v/po > 0. 

(We omit the details of the calculation, since the procedure is the same as in 
[^] and is straightforward but so tedious.) This confirms the existence of the 
global solution. □ 

The solution (p) can be caluculated from (|3.38|) , and in the region |T| < s'^Z, 
it takes the constant value, 

(p) = po f 1 + 77^1 >Po. (4.29) 
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A,5 = A,6 

A,2 = A,3 ■ 







Uo-lVpo 



-Uo-lVpo 



Uo+lVpo 



-Uo+lVpo 



T 

Xi = Xi 



Figure 14: Evolution of the Riemann invariants for < uq < 



This explains the elevation of the pulse level observed in Fig.^. In Fig.|T4|, 
we plot the evolution of the Riemann invariants. The corresponding solution 
p{T, Z) at Z = 5(3 ^ 1.58 of the numerical simulation of the NLS equation 
( pTTD is shown in Fig.||, where po = l,uo = 3(3 ^ 0.95 with (32 = —0.1. Note 
in Fig.^ that the non-oscillating part with the constant level ( [4.29 ) in the 
center region |T| < s'^Z corresponds to the g = 2 (degenerated into g = 0) 
solution, and the oscillating parts in the regions s'^Z < \T\ < Z to the 



ones (degenertated into (? = 1), as predicted in Fig.^ One should compare 
Fig.^ with Fig.|13| to see the difference in the center regions. 



(Hi) The case with > uq > 

In this case, the initial data gives decreasing functions for both Ai and A2 
(see Fig.^D, and we have: 

Theorem 4.3 For > uq > —2J^, the g = NLS-Whitham equation 



admits the global solution with the initial data M-lj ). 
Proof. In this case, Ai and A2 can be regularized as 



Ai(T) = lim Ai(r 



e), and A2(r) 



lim A2(T 



(4.30) 
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Figure 15: Evolution of the Riemann invariants for > Uq > — 



Then from Eq. (|4.2| ) of Sk for g = 0, we obtain 



-Si' 



-UO + ^/pO > 0, 



Thus we have S2 > S2 > > Si > , which imphes the assertion of the 
theorem. □ 

We also find the solutions Ai, A2 and (p(T, Z)) = p(T, Z): For T > 0, 

for T > S2Z, 



(4.32) 



uo - 2v/po, 

Xi{T,Z) = \ -±T-^{uo + 2^o), for s+ZkTks^Z, 
-Mo - 2^/p^, for < T < S2Z, 



X2{T,Z) = uo + 2^, for VT > 0, 



and 



PiT, Z) 



Po, 



for T > Sr, Z 



2 ^) 



9^2 



{T+(no + 2Vp^)Z}^ for 4Z < T < ^2 ^, 



(4.33) 



Po 1 + 



Mo 



2Vp^y ' 



for < T < 4Z, 



32 



2 




Figure 16: Deformation of the pulse p(T, Z) for > Mq > —l^fp^: Solid 
curve, the numerical result; Dashed curve, the solution of Eq.( [4.33| ). 



For T < 0, note the symmetry Ai(— T, Z) = —\2(T,Z) and p{—T,Z) = 
p{T,Z). These solutions are shown in Fig.|15] and Fig.p!6|. In Fig.|16|, we also 
plot the numerical result p{T,Z) at Z = 10/9 ~ 3.16 of the NLS equation 
(|2.1| ) where po = 1,uq = —3/3 ~ — 0.95,/32 = —0.1, which agrees quite well 
with the analytical solution ( [4.33[ ). Note that the level of p for |T| < s^Z has 
the same formula as the previous case except the sign of Uq. This explains 



the bit-wise depression in the pulse level shown in Fig.|Tl|, and it leads to the 
deformation into a RZ-like pulse observed in the experiment 0. This also 
gives a theoretical limit of the strength of the initial chirp, i.e. uq > —2y/po, 
to avoid a destruction of the pulse. In fact, we will see below that the case 
with Mo < —l^fp^ the level of the signal becomes zero starting from T = 0, 
the point of the discontinuity in chirp. 

(iv) The case where Uq < —2yfp^: 
We first have: 

Lemma 4.4 The initial data ^■'^ ) for Uq < — 2y/p^ can be considered as a 
g = I initial data where the Riemann invariants Ai, ■ ■ ■ , A4 have the following 
values, 



A = / "0 + '^y/Po, for T <0, 
' \uo- 2^, for T>0, 
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Figure 17: Evolution of the Riemann invariants for uq < —IJ^ 



A2 = Mo + 2yp^, for VT, 



(4.34) 



A3 = -Mo - 2v^, for VT, 



A4 



-Mo + 2^, jor T <0, 
-uo-2^, for T>0. 



Proof. The calculation to show (p) = p and (pn) = pu is the same as the 
previous cases, and we omit it. □ 

In Fig.0, we plot the initial data as the solid lines. Since all the Riemann 
invariants Afc satisfy the conditions in Corollary ^4.1|, we obtain: 



Theorem 4.4 For uq < — 2y/p^, the g = 1 NLS-Whitham equation with the 
initial data (|^. has a global solution. 



Proof. We determine the evolutions for Ai and A4. The regularization of 
Ai and A4 at T = is given by imposing 



Ai(T) = ^hm Ai(T - e), and A4(T) = ^lim A4(T - e). (4.35) 
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Figure 18: Deformation of the pulse p{T, Z) for uq < —2J^ 



Then we have > , and therefore we obtain the global solution as in the 
previous cases. In fact, from the formula of s^, we have 



-sX = -Uq + ^ > 0, 

(4.36) 

-s" = -uq - 2y/p^ > 0. 



This completes the proof. □ 

We note in particular that the solution (p(T, Z)) in the center region |T| < 
S4Z becomes zero, which is obtained from Eq. (|3.38|) with ai = 0, (T2 = 
-2X1 ^1 = ^3^1 and ai = -A^. In Fig.|T8|, we plot the pulse profile p{T, Z) 



at Z = 3/5 ~ 0.95 obtained by the numerical simulation of the NLS equation 
( P?TD with po = Ij'Wo = —9/3 ^ —2.85 and (32 = —0.1. As was predicted by 
the theory, the pulse level at T = becomes zero right after the propagation. 
One should also note that the genus is actually degenerated into = for 
VT, even though we need the g = 1 regularized initial data for the global 
solution. This implies that the solution is indeed g = 1, but there is no 
oscillations in the solution. The small oscillations observed in Fig.|l^ are due 
to a nonzero value of (32-, and disappear in the limit /?2 — »■ 0. 

In this Section, we have considered the evolution of the Riemann invari- 



35 



Z < 



z = o 



z >0 




Figure 19: Deformation of the Riemann surface for < uq < 



ants for Z > in the cases with several different values of initial chirp uq. 
We then recognized that some of the cases gave degenerate genus initial data, 
and choosing appropriate genus, we found the global solutions for the corre- 
sponding NLS-Whitham equations. We also note that the original system of 
equations (|2.2|) and (|2.3| ) is invariant under the change of variables Z —Z 
and u — >■ —u. Hence for example the evolution of the case [ii) for Z > 
can be extended to the evolution for Z < with the case [in) under the 
change Z —Z and Uq — * —Uq. Thus we have described a change of genus 
along the Z-axis, and in the example we see the Riemann surface of g = 
for Z < and of g = 2 (degenerated into g = 1) for Z > 0. This process is 



illustrated in Fig. 19 



5 Effects of the higher order dispersion 

As shown in |TT|, |T^, the important higher order effects in a long-distance 
and high-speed optical communication system are described by the perturbed 
NLS equation, 

+ y ^ + 1^1 ? = ^ + ^7^(kl q) + -nq^\q\ , (5-1) 

where jS^, 7, are real constants describing the third order linear dispersion, 
nonlinear dispersion due to the delay effect of the Kerr nonlinearity and the 
Raman (dissipative) effect. Since the NRZ system operates in a small GVD 
{P2) regime, the third order dispersion (/^s) becomes particularly important. 
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Also the Raman term may not be so important in the regime of pulsewidth 
with scores of pico-second. So we consider the effects of /?3 and 7 in Eq. (^.l|) . 
In a real system, we also use a region where /^s becomes (or make) so small, 
and assume that /^s oc (/32)^''^- Then as shown in Eg. (^.1|) with = 
can be Lie-transformed to an integrable system (the NLS + its hierarchy), 

'dZ ^ Ydf^ + \l\l = 'Ydf^ df ^ ^ orders). (5.2) 

In a small (32, the transformation is however valid only for a weak nonlinear 
case, e.g. |gp ~ 0{P2)- Equation (|5.2|) without the higher orders has also 
infinite number of conservation laws with the same conserved densities as the 
case of the NLS equation. In fact, this can be also written in the Lax form, 

(9/ 

-t— = [B,L], with B = B2+(3B,, (5.3) 

where S„ are defined in (U), and /3 = -13^/{2P^) with /^s = The 
explicit formula of B^ is given by 

B;=^-D'' + pD + pu + ]^Dp. (5.4) 

Now let us take the limit /5 — >• 0, and assume both p and u are smooth. 
Then as in the previous section we have the g = Q approximation of the equa- 
tion (|5.2|) , which is obtained in the form similar to ( 3.101) with the rescaling 

^ = ^(Q2 + /3Q3), (5.5) 

where Q3 = [A^]>o = P^/3 + pP + pu. In terms of p and m, Eq.(p.5|) gives a 
system. 



^ ( P\ _ ( u + (3{u^ + 2p) p + (3{2up) \ 9 fp 



dZ\u) \ 1 + /?(2m) u + p{u'' + 2p)) dT 



u 



(5.6) 



Since the Lax operator for the equation ( p.2|) is the same as the NLS case, 
the spectral curve in this limit is given by (|3.8|) , and the system ( |5.6| ) is also 
referred as the (7 = Whitham equation. The eigenvalues of the coefficient 
matrix of (|5.6| ) are then given by 

si,2 = n + /3(M2 + 2p) T(l + 2/3u)v/p (5.7) 
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which give the characteristic velocities of the Riemann invariants Ai and A2 
in (I2.10I ), i.e. Ai,2 = m =F 2^^. In terms of the Riemann invariants, Eqs. (|5.7|) 
are expressed as 

s, = ^(3Ai + A2) + |(5A2 + A^ + 2A1A2), 

S2 = ^(Ai + 3A2) + |(A? + 5\l + 2A1A2), 

Because of the nonhnear terms of Ai and A2, the characteristic velocity does 
not satisfy the double monotonicity in Corollary [4.1| . This implies that the 
g = Whitham equation ( ^.6|) develops a shock even in the case where 
the Riemann invariants are both decreasing functions. To demonstrate the 
situation, we study the initial value problem of (|5.6|) with the initial data 
given by (|2.8|) . As in the previous section, we calculate the characteristic 
velocities of the pulse edges at T = — Tq and at T = Tq (see Fig.^). 

Let us first compute the velocity near T = Tq, where the Riemann invari- 
ant A2 decreases from 2y/p^ to — 2y/po, and Ai = — 2^/po- From Eq.( |5.8|) , S2 
near T = Tq is given by 

S2{Xi = -2y/p-o, A2) = ^A^ + ^(3 - 2/3^)A2 - - ^Vp~o). (5.9) 

To avoid a shock, we require that S2 is a monotone increasing function in 
A2 G [—2y/po, 2y/po]. This implies that for /3 > the velocity S2 has to take 
the minimum in the region A2 < —2J^, and we have the relation. 



=-3/3v/p^+^>0. (5.10) 



dX2 

Similarly, for P < 0, S2 has to take the maximum in A2 > 2^/po, and 

l?L ,^ = ^^Vp'o + 1>0. (5.11) 
From these results we obtain 

<f3<-^. (5.12) 



Po 4^po 
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To compute also the velocity si near T = — Tq, we note the symmetry of the 
Riemann invariant form of the g = Whitham equation (|5.6|) , that is, for 
the reflection T —T the form is invariant under the change of variables 
Afc — > —X-s-k for k = 1,2 and f3 — > —/3. Therefore we obtain the condition to 
avoid a shock near T = —Tq, 

^ <P<^—. (5.13) 



Thus we have : 



Proposition 5.1 The system ( \6.(^ with the initial data l \2.^ admits a global 
solution if the following condition is satisfied, 



This provides a condition on design parameters in the second and third order 
dispersions and the power level for a smooth NRZ pulse propagation. 

The solution of ( |5.6| ) can be obtained as follows: Following the proof of 



Proposition |2.1| , for example near T = Tq, we have a single equation for A2, 

S2i-2Vp-o,X.)^, (5.15) 



dZ ^ ' dT 

where S2(— 2y^, A2) is given by (|5.9|) . The solution of this equation can be 
obtained in the hodograph form, 

T + S2i-2^oA2)Z = To, (5.16) 

from which A2 gives a monotone decreasing function in T from A2 = 2y/p^ to 
A2 = —2y/p^, if the condition ( |5.14| ) is satisfied. In Fig.^, we plot the pulse 



profiles piT,Z) at Z = 5/3 ~ 1.58 obtained by the numerical simulation 
of (|5.2|), where we used (32 = —0.1 and = 0.02, giving a case of shock 
formation i.e. (3 = 0.63 > 0.25 with po = 1- Iii this example, we predict the 
shock formation from both edges T = ±To, which are due to the violation of 
the conditions (|5.12| ) and ( p.l3| ). The asymmetric deformation of the pulse 



due to the third order dispersion has been observed in the experiment 



It is also interesting to study how one can regularize where the 



condition (|5.14|) is violated. For example, if /5 = 3/(2i/po), then S2 takes 
minimum at A2 = 0. So A2 becomes multi-valued function in T for Z > 0, 
which indicates a shock formation, and one needs to regularize the initial data 
by choosing an appropriate genus. This will be further discussed elsewhere. 
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Figure 20: Pulse profile under the effect of third order dispersion. 



6 WDM-NRZ system 

In a WDM system, signal in each channel has a different carrier frequency, 
and the model equation of the system can be obtained by setting the electric 
field q of the NLS equation in the form. 



N 



(6.1) 



Here qj represents the electric field in the j-th channel having carrier fre- 
quency, say Vtj. Because of the difference in the carrier frequencies, pulse in 
one channel has different group velocity from those of the others, and we see 
pulse-pulse collision between channels. As a simple but fundamental case, 
we consider here the 2-channel WDM system. Then substitutng ( |6.1|) into 
the NLS equation and ignoring the mismatch of the frequency in the 



FWM terms, we obtain the following coupled NLS equation |21 



.dqi 
.dq2 



f3'^ d'^qi , /, |2 I I |2^ 



(6.2) 



'dZ 2 dT^ 



+ (l'?2p + (y\qi\^)q2 
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where a = 2. The vahdity of this model has been discussed in where a 
dispersion management technology is used to reduce the FWM and to keep a 
small GVD. We should also like to mention that the system (|6.2D with a = 1 



is integrable, and is called the Manakov equation [20 



Let us first summarize the result in 17 for the 2-channel WDM case 



[a 



2). The main result in ||T^ is that there exists a critical frequency separation 
AQc such that the system with the channel separation 1^21 — ^21 < 
leads to a hydrodynamic-type instability. It is also interesting to note that 
this instability appears even in the case of a = 1, and is not essential for 
the specific value of a. So we study this instability based on the Manakov 
equation, which enables us to characterize the instability as a deformation of 
the branch points on the Riemann surface defined by the spectral curve of 
the Lax operator for the equation. 

The Lax operator Lm for the Manakov equation takes a similar form as 
of the NLS equation, (pTip, and is given by 



D + 



where pk and Uk are of course defined by 



qk{T, Z) 



PkiT, Z) exp U-^(Tfc(T, Z) 



for k = l.2. 



(6.3) 



(6.4) 



As seen from this formula, the Manakov equation can be easily extended to 
an integrable iV-coupled NLS equation which may correspond to a model 
equation of A^-channel WDM system. Then the Manakov equation can be 
written in the Lax form (p.5[) with B2 = {L\.j)+/2, i.e. 



dL 



M 



dZ 



[B2,L 



M\ 



(6.5) 



where we have rescaled Z to Z/jS. 

Now taking the limit (3 —>■ and ignoring all the terms including (3, we 
obtain a hydrodynamic-type equation. 



(6.6) 





(Pi] 




1 Ui 


Pi 





^) 




(Pi\ 


d 


Ui 




1 


Ml 


a 





d 


Ui 


dZ 


P2 












P2 


df 


P2 








\ a 





1 


U2) 
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Note that this equation can also be derived from the original system ( |6.2| ) 
with arbitrary a by substituting ( |6.4|) into (|6.2|) in the same approximation. 
However a system with a ^ 1 may not be integrable and the corresponding 
Whitham equation may not be constructed in a systematic manner. 

To see the hyperbolicity of ( |6.6|) , we first calculate the characteristic poly- 
nomial of the coefficient matrix in (|6.6|) , 

[{u, - sf - pi] [{U2 - sf - P2] - aViP2 = 0. (6.7) 

As a practical example, we set the equal intensity p\ = p2 = po and the 
opposite frequency shifts Ui = —U2 = Uq, then we have a condition for the 
hyperbolicity of (|6.6|) , that is, all the eigenvalues of the coefficient matrix are 
real, if the following condition is satisfied, 

\uo\ > [(l + a)po]'/'. (6.8) 

This then gives a minimum frequency separation for 2-channel WDM sys- 
tem for a stable NRZ pulse propagation. The instability corresponding to 
the complex eigenvalues (characteristic velocity) is called hydrodynamic-type 
instability in [jl^. Thus the instability is due to the change of type of the 



system ( |6.6|) from hyperbolic to elliptic. 

We now consider the geometry of the instability using the Manakov equa- 
tion, and discuss the associated Whitham equation to describe a schock for- 
mation in WDM system. In the limit (3 0, we again assume the wave 
function of Lm4' = -^0 iii the WKB form ( p.7|) , then the spectral curve 
A = A(P) is given by 

A = P + 7^ + 7J^, (6.9) 

where P = dS/dT, and the evolution of P is given by dP/dZ = dQ2/dT 
with Q2 = (l/2)[A2]>o = PV2+P1+P2, the limit of ^20/0 in (§3p. In Fig.^, 



we plot the graph of Eq. ( |6.9|) on the P - X plane. As we see from Fig.^ there 
are at most three roots of P in Eq.( |6.9| ) for each A, defining the three sheeted 
(non-hyperelliptic) Riemann surface. The characteristic polynomial (|6.71) is 
then given by the equation dX/dP = 0, and the eigenvalues are obtained 
from ( |6.9| ) at the roots of this equation. Thus the spectral of Lm of ( |6.3| ) is 
given by 

Spec(LAf) = (-00, Ai] U [A2, A3] U [A4, 00). (6.10) 
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Figure 21: The g — algebraic curve (6.9). 




Figure 22: The g — Riemann surface corresponding to the curve (6.9). 
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We also note that the spectral curve ( |6.9| ) defines the g = Riemann surface 
whose topology is a sphere obtained by a compactified P-plane (see Fig .[2^). 

To see these more precisely, let us study a detail of the algebraic curve 
(|6.9|) . We first note Eq.(|6.9|) can be written as 



V-3«o(A)i/-/9o(A) = 0, 



(6.11) 



where 2y := 3P — {\ + Ui + U2), and ao,Po are the monic polynomials of 
degree 2 and 3 in A, 



"0 



- {ui + U2)X + ul + ul- U1U2 - 3(pi + P2) 



(6.12) 



9 

--(piMl + P2U2) + 9(piM2 + P2U1). 



+ u\- AU1U2 + 3(pi + P2)} A 



The roots of ( |6.11 ) are given by 



for 



0,1,2, 



(6.13) 



(6.14) 



where a; = exp(27ri/3), and are given by 



(6.15) 



Each root lives in a different sheet, and in particular the root marked by 
£ = in ( |6.14 ) lives in the fundamental sheet having the property P ^ A as 
A ^ 00. With this root. Eg. ( |6.6|) can be identified as the g = Whitham 
equation in the form ( |3.21| ) where the differentials uji and 002 are given by 



u;, = dP = Ul 



1 



+ 



1 



3(z/+)2/3 d\ 3(z/_)2/3 d\ 

dX, 



dX 



1 + ^ + 



A2 



UJ2 = dQ2 = PdP 



MlPl + U2P2 , 

A+ + 



dX. 



(6.16) 
(6.17) 
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The branch points (Riemann invariants) A^, k = 1, ■ ■ ■ ,4 are obtained from 
the descriminant Ao(A) of the equation (|6.11|) , that is, 



4 



{ui - 1[{X - Xk) = 0. (6.18) 



k=l 



X 



2 



In the case where ui = —U2 = uq and pi = p2 = po, the roots of Eq. (|6.18| ) 
are given by 

\ {24 + IOuIpo -pI± (Aul + Po)'/'v^} . (6.19) 

Uq 

We also culculate the characteristic velocities Sk = Pk = Pi^k) which is given 
by the roots of dX/dP = 0, 

P' = ul + po± ^{Aul + po)po. (6.20) 

If Uq > 2po of (|678|), these roots are all real. Note that = A(Pfc) with the 

order ings Ai < A2 < A3 < A4 and Pi < P2 < -P3 < -Ri- Then the condition 
( |6.8D indicates the non-closing of the interval [A2, A3]. 

Now let us consider the case with nonzero genus. The corresponding 
algebraic curve may be given by 

Ay' - 3ahiX)y - f3hiX) = 0, (6.21) 

with 



2h+2 

ah= n (-^ ~ ^i)' 

3/1+3 

Ph= Ui^-Vk), 

k=l 



(6.22) 



where we have the relation 3X]j=t^0 ~ '^Y.t^^ Vk- Then the descrimenant 
Ah := — (4/27)(/3|— a|) of the curve ( |6.21D gives a polynomial of degree Gh+A 
in A, and it indicates that there are at most 3h numbers of genus openning 
in the Whitham equation. Note that the genus of the (non-hyperelliptic) 
Riemann surface defined by the curve ( |6.21| ) is generically given by 

g = 3h. (6.23) 
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The number 3 implies that gap-opening appears in between 2 sheets of the 
total 3 sheets. A further discussion on this subject will be given elsewhere. 
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Figure Captions 

1. Various types of a 16 bits-coded signal (0010110010111100). 

2. Evolution of the Riemann invariants. Doted curve are the initial data; 
Solid curves are the solutions Xk(T, Z) for Z > 0. 

3. Deformation of NRZ pulse (Dam breaking problem). The dotted curve 
in p shows the numerical result of the NLS equation. 

4. Optical shock due to initial chirp. 

5. The (? = algebraic curve (|3.8| ). 

6. The g = Q Riemann surface corresponding to the curve ( p.8| ). 

7. The g = 1 algebraic curve y"^ = Rg{\) in ( |3.26|) . 

8. The 2-sheeted Riemann surface with branch cuts and canonical cycles. 

9. The function ^^(A) with i = 2k + 1. 

10. The initial phase modulation for one-bit pulse. 

11. Evolution of NRZ pulse with initial phase modulation, a) Without, 
and b) with the initial modulation. 

12. Evolution of the Riemann invariants for mq > 2y/pJ^. 
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13. Deformation of the pulse p(T, Z) for uq > 2y/po- 

14. Evolution of the Riemann invariants for < Mq < 2y/p^. 

15. Evolution of the Riemann invariants for > mq > —IJ^. 



16. Deformation of the pulse p(T, Z) for > Mq > —2^fp^. Solid curve, the 
numerical result; Dashed curve, the solution (4.34). 



17. Evolution of the Riemann invariants for mq < — 2y/p^. 

18. Deformation of the pulse p(T, Z) for uq < — 

19. Deformation of the Riemann surface for < Mq < '^^[Pq- 

20. Pulse profile under the effect of third order dispersion. 

21. The (7 = algebraic curve (|6.9|). 



22. The (7 = Riemann suface corresponding to the curve 
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